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ABSTRACT 

Anomalous microwave emission at 20-40 GHz has been detected across our Galactic 
sky. It is highly correlated with thermal dust emission and hence it is thought to be 
due to spinning dust grains. Alternatively, this emission could be due to synchrotron 
radiation with a flattening (hard) spectral index. We cross-correlate synchrotron, free- 
free and thermal dust templates with the WMAP 7-year maps using synchrotron 
templates at both 408 MHz and 2.3 GHz to assess the amount of flat synchrotron 
emission that is present, and the impact that this has on the correlations with the 
other components. We find that there is only a small amount of flattening visible 
in the synchrotron spectral indices by 2.3 GHz, of around A/3 0.05, and that the 
significant level of dust-correlated emission in the lowest WMAP bands is largely 
unaffected by the choice of synchrotron template, particularly at high latitudes (it 
decreases by only ~7 per cent when using 2.3 GHz rather than 408 MHz). This agrees 
with expectation if the bulk of the anomalous emission is generated by spinning dust 
grains. 

Key words: diffuse radiation - Galaxy: general - radiation mechanisms: general - 
radio continuum: ISM 



1 INTRODUCTION 

Maps of high-frequency (> 10 GHz) radio emission contain 
multiple distinct astrophysical components that must be dis- 
entangled to study both the astrophysics of our Galaxy and 
the cosmological information contained in the Cosmic Mi- 
crowave Background (CMB) radiation. In particular, multi- 
frequency microwave instruments such as the Wilkinson Mi- 



crowave Anisotropy Probe ( WMAP, Bennett et al. 2003a I 
and the Planck satellite ( Tauber et al.|2010 ') primarily mea- 
sure combinations of synchrotron emission from electrons 
spiralling around magnetic fields (low frequencies); free- free 
emission from the interaction of electrons with ions (interme- 
diate frequencies); anomalous microwave emission (AME) 
thought to be from spinning dust grains (peaking at ~20- 



40 GHz; see e.g. Draine & Lazarian||1998, ,de Oliveira-Costa 


et al.|2004 


IWatson et al.|2005| Ysard, Miville-Deschenes & 


Verstraete 


2010 


Planck Collaboration et al. 


2011 


I; CMB 



emission (intermediate to high frequencies); and thermal 
dust emission (high frequencies). 

These various components can be separated from each 
other by fitting templates of the different emission mech- 
anisms to the data, i.e. using maps of single components 



to separate the high frequency maps based on the emis- 
sion morphology (see e.g. [Kogut et al.||1996{ [Banday et al.| 



2003 



Davies et al. [200^ 1. However, problems arise when the 
morphology of the emission changes with frequency, which 
can result from the differing spectral indices of the various 
astrophysical objects that are generating that emission. In 
particular, there is a wide range of observed synchrotron 
spectral indices, and the overall behaviour of the emission 
may change at higher frequencies - either steepening due to 
synchrotron ageing, or flattening due to a recent injection of 
energetic electrons or multiple components of emission. 



In this context, [Bennett et al.| ( |2003b[ ) have suggested 
that the dust-correlated emission observed by WMAP could 
be due to flat (hard) synchrotron radiation (with a spectral 
index /3 ~ —2.5) that is not well traced by the standard 
low-frequency map of synchrotron emission, as opposed to 
being due to a different emission mechanism such as spin- 
ning dust. These two scenarios are sketched out in Figure [l] 
In this paper, we apply the template fltting method to the 
WMAP 7-year maps using synchrotron templates at both 
408 MHz and 2.3 GHz. The 408 MHz map is the standard 
tracer of synchrotron radiation used in this type of analy- 
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sis, however due to its low frequency it will predominantly 
trace steep spectrum radiation, whereas the 2.3 GHz map at 
nearly 6 times higher frequency will trace more of flatter- 
spectrum radiation. We look for an increased correlation 
between the WMAP maps and the higher frequency syn- 
chrotron template, and a corresponding reduced correlation 
with the dust template, which would be indicative of the 
anomalous emission being flat spectrum synchrotron radia- 
tion. 

In order to judge how sensitive the 2.3 GHz map might 
be to flat-spectrum synchrotron compared with the 408 MHz 
map, we use the r.m.s. values of the 408 MHz and FDS maps 
(5.9 K and 6.8 (xK respectively, see |Davies ct al. 2006) com- 
bined with the template coefficients measured later in this 
paper, and assume spectral indices of /3 = —3.01 for steep- 
spectrum emission and —2.5 for the flat spectrum emission. 
Assuming that all of the dust-correlated emission measured 
at 22.8 GHz is flat spectrum emission, and assuming no spec- 
tral curvature in either the steep- or flat-spectrum compo- 
nent, then it would constitute 29 per cent of the 2.3 GHz 
emission, compared with 14 per cent of the 0.408 GHz emis- 
sion. Thus, we would expect to see a considerable difference 
in the template coefficients. Given this, one might expect 
the 22.8 GHz dust template coefficient to change by a factor 
~ 14/29 if the dust-correlated component were completely 
due to a flatter spectrum synchrotron spectrum. Though 
there is evidence of curvature in the synchrotron spectrum, 
these flgures give a qualitative idea of the difference in rela- 
tive strengths of the components in the low-frequency tem- 
plates that motivates this work. 

In Section [2] we discuss the data used here, and the 
choice of the specific foreground templates. We briefly sum- 
marise the cross-correlation method that we use to carry out 
the template fitting in Section [3] Our results are presented 
in Section [4] along with an assessment of their robustness 
as a function of resolution. Galactic masks, and different 
template options. This section also looks at a few selected 
areas of Galactic emission, where we confirm that our results 
are valid both for the whole sky and for sections of it. We 
conclude in Section (5] 



2 DATA AND TEMPLATES 

For this analysis, we use a combination of low frequency 
maps from ground-based radio telescopes which map the 
synchrotron radiation (described in section 2.1); optical 
Hq emission maps tracing the free-free emission (section 
2.2) and frequency-scaled maps of the dust emission from 
IRAS (section 2.3) as our foreground templates, and which 
we then fit to the maps from WMAP (section 2.4). The 
properties of the maps used here are summarised in table [l] 
All of the maps are in the HEALPl>Qpixelization scheme 
( Gorski et al.|2005| >. The maps are smoothed to the desired 
resolution (nominally 3°) whilst in their original A^sido (ei- 
ther 512 or 256) before being down-sampled to a nominal 
Aside = 64 (i.e. a pixel size of 55 arcmin). We mask all 
of the maps using the extent of the 2.3 GHz observations. 
In addition, we use a combination of WMAP Galactic plane 
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Figure 1. A sketch indicating the different emission mechanisms 
present at various frequencies. Spectral shapes are plotted in ar- 
bitrary flux density units, with spectral indices a = 13 + 2. Two 
possible synchrotron components are shown in black: one steep 
spectrum (S oc i^"^-"), the other flat spectrum but turning over 
at WMAP frequencies (S oc /[l + [u/A5GHz]^]). An illustra- 

tive spinning dust model is also shown in grey. Free-free emission 
{S oc is shown in red, and thermal dust {S oc u^'^) in 

magenta. The cause of the anomalous emission at WMAP fre- 
quencies could be due to either spinning dust or flat synchrotron 
emission that turns over at WMAP frequencies. 



and point source masks to exclude strong emission that may 
bias the results (particularly due to the absorption of Hq by 
thermal dust). 

Throughout this work we use the spectral index con- 
vention of r oc . We use Rayleigh- Jeans (RJ) antenna 
temperatures throughout. Where maps are in CMB tem- 
peratures we convert to antenna temperatures using Trj = 
TcMBX^e^/{e^ — 1)^, where x — hu/kTcMB, in which 



TcMB = 2.725 K. 



2.1 Synchrotron templates 

The standard all-sky synchrotron map, which is used 
in a wide variety of component separation analyses, is 
the 408 MHz map (IHaslam et al.l 119811 119821), hereafter 



sync[H0.4]. This is currently the only high resolution (~ 1°) 
map of the entire sky that is available at a frequency of less 
than a gigahertz. It combines observations made with four 
large single-dish telescopes (Jodrell Bank Mkl and MklA, 
Effelsberg and Parkes) and is absolutely calibrated using the 
404 MHz survey by |Pauliny-Toth fc Shakeshaftl |T962| (ex- 



trapolated between 404 and 408 MHz using a spectral index 
of —2.5). The 408 MHz map hence has an accuracy of ~ 3K 
in absolute temperature and ±10 per cent in Tb scale. 

There are a number of different versions of this avail- 
able (see e.g. |Davies, Watson & Gutierrez| |1996| [Plata 



nia et 



"ar|[2003f 



We use the version available from the 
Legacy Archive for Microwave Background Data Analysis 
(LAMBDAQ This version has been de-striped and point 
sources have been removed. Although the resolution of the 



http : //healpix • jpl • nasa . gov 



http : //lambda . gsf c. nasa . gov/product /f oreground/haslam_ 



|408 . cf m| 
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Table 1. The properties of the template maps (top) and input 
maps (bottom) used in this analysis. Where maps have been pre- 
smoothed, we provide both the smoothed and raw resolutions. 
Hcf is used as a proxy for free-free emission, so we do not quote a 
frequency for it here. The thermal noise per observation, cr, is rele- 
vant for the WMAP maps; the others are systematics-dominated. 



Name 


Freq. (GHz) 


Res. 


CTO (RJ) 


Notes 


sync [HO. 4] 


0.408 


60' 




orig. 51' 


sync[J2.3] 


2.326 


60' 




orig. 21' 


dust[FDS94] 


94 


6' 






Ha[F03] 




-6' 




Mixed res. 


Ha[DDD] 




60' 




/d = 


Ha[DDD0.33] 




60' 




fd = 0.33 


WMAP K 


22.8 


60' 


1.418 niK 


orig. 0?88 


WMAP Ka 


33.0 


60' 


1.429 mK 


orig. 0?66 


WMAP Q 


40.7 


60' 


2.105 mK 


orig. 0?51 


WMAP V 


60.8 


60' 


2.854 mK 


orig. 0?35 


WMAP W 


93.5 


60' 


5.253 mK 


orig. 0?22 



survey is nominally 0?85, the post-processing that has been 
applied to it has smoothed the map slightly; we assume that 
the map now has 1° resolution. As we are looking at the 
diffuse emission, which is not significantly affected by the 
precise beam size, this assumption is sufficient for our anal- 
ysis. 

At 5.7 times higher frequency than sync 



H0.4], the 

2.326 GHz map of Jonas, Baart & Nicolson ( 1998 1, hereafter 



sync[J2.3], made from observations with the 26-m HartRAO 
telescope, can also be used as a synchrotron template. Al- 
though not a full sky map, this still covers 67 per cent of the 
sky at high resolution (20 arcmin) and high sensitivity. The 
map was made differentially, with an arbitrary zero point 
that was determined using the 2.0 GHz absolute measure- 



ments by Bersanelli et al. ( 1994 1 (assuming a spectral index 
between 2.0 and 2.3 GHz of —2.75), and an uncertainty in Tb 
of ~5 per cent. As such, the map should be completely inde- 
pendent of sync[H0.4], and in the presence of flat-spectrum 
synchrotron emission it should be a better tracer of the high- 
frequency synchrotron that is seen by WMAP. 

The two synchrotron maps are shown in Figure |2] The 
change in the morphology of the synchrotron emission be- 
tween these two frequencies is obvious. For example, the 
north polar spur has a steep synchrotron spectrum and as 
such it is much less noticeable in sync[J2.3] than sync[H0.4]. 
We note that a 1.4 GHz map also exists ( Reich|[l982[ 



Reich fc Reich|1986| [Reich, Testori fc Reich|2001[ ), however 
we do not make use of this here as the higher frequency map 
will provide a better tracer of flat spectrum synchrotron 
radiation. In the future, this work can be repeated using 
a 5 GHz map from C-BASS (|King et al.|[2010l), which will 



provide an even better template of flat spectrum synchrotron 
radiation. 



2.2 Free- free templates 

We use a map of the Ha line emission at 656.28 nm in units 
of Reyleigh (R) as a proxy for the free-free emission, as is 
conventional for component separation due to the difficulty 
in cleanly separating large-scale free-free from synchrotron 
emission at radio frequencies. There are, however, a number 



of issues with any map of Ha emission - including geocoro- 
nal emission from the Earth's upper atmosphere, baseline 
effects, stellar residuals and dust absorption. There are two 
full-sky Hq maps available, both of which attempt to mit- 
igate these effects, namely those by [Dickinson, Davies fc| 



Davis (20031 and Finkbeiner (20031. Both of these are com- 



posites of several different large-scale ground-based surveys. 

The map produced by 'Finkbeiner! ('2003'), hereafter 
Hq:[F03], combines data from the Wisconsin H-alpha Map- 
per survey (WHAM; [Reynolds et al.|1998|l, the Southe rn H 
Alpha Sky Survey Atlas (SHASSA;|Gaustad et al.[2(30T|) and 



Virginia Tech Spectral line Survey (VTSS; Dennison, Simon- 
|etti fc To pasna 1998). Unfortunately, these three surveys 
have different resolutions and fields of view, which means 
that the combined map cannot be easily smoothed to a single 
resolution. We assume that the resulting map has a resolu- 
tion of 6 arcmin, which will lead to over-smoothing in parts 
of the map where the observations were made at a lower res- 
olution (see Ghosh et al.pOll for a power-spectrum analysis 
of this issue). Again, as we are focusing on the diffuse emis- 
sion, this should not have a significant effect on our results. 
Additionally, most of the fields where this is the case are 
in the Northern hemisphere; of the Southern area that we 
are using, only 10 per cent of the sync[J2.3] survey area is 
affected (or 7 per cent when the KQ85 mask is applied). 
The maps produced by [Dickinson, Davies fc Davis| 



( 2003 1 avoid this resolution issue by combining only WHAM 
and SHASSA with suitable smoothing to a resolution of 1°. 
We use two versions of this map - one map that has been 
corrected for dust absorption of the Ha emission using the 



100 \im dust map of Schlegel, Finkbeiner fc Davis ( 1998 1 



with an effective fraction of the dust that contributes to- 
wards Ha absorption of /d = 0.33 (Ha[DDD0.33]), and a 
second map without this correction factor (Ha[DDD]) for 
more direct comparison with the results using Ha[F03]. 

We use the Ha[DDD] map as our standard map, 
and will investigate the effects of using the Ha[F03] and 
Ha[DDD0.33] maps in section 4.2. We correct the syn- 
chrotron maps for free-free contamination using the 
Ha maps. We use the values presented in [Dickinson, Davies[ 



fc Davis ( 2003 1 assuming that the gas has an electron tem- 
perature Te = 7000 K everywhere on the sky (although 
the electron temperature is known to vary across the sky, 
this has yet to be quantified and mapped). At 408 MHz 
this correction factor is 51.2 mK/R and at 2.326 GHz it is 
1.33 mK/R. The inclusion of these free-free corrections to 
the synchrotron templates does not significantly alter our 
later results. 



2.3 Thermal dust 



We use the thermal dust map created by [Finkbeiner, [ 
[Davis fc Schlegel[ ( [19991 hereafter dust[FDS94]) via an ex- 
trapolation of the 100 ^m IRAS and 100/240 ^m COBE- 
DIRBE maps to 94 GHz. This map has a resolution of 6 ar- 
cmin. We use the version available from LAMBD^ In the 
future, maps of the cold thermal dust as seen by Planck HFI 



http : //lambda . gsf c. nasa . gov/product /foreground/ dust. 



|map . cf m[ 
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-O.B0 ^^^^^^K i^^^H 1.? Log (R) 0.10 ^^^^^H ^^^H S.l Log (uK) 

Figure 2. The input templates, smoothed to 3°, using A^sidc = 64 and plotted with a logarithmic colour scheme: sync[H0.4] (top-left), 
sync[J2.3] (top-right), Ha[DDD] (bottom-left) and dust[FDS94] (bottom-right). The maps are in MoUweide projection and are in Galactic 
coordinates with longitude increasing to the left. All maps have been masked with the 2.3 GHz survey area and the KQ85 mask. The 
morphological differences between the different emission mechanisms can be clearly seen. 



at 857 or 545 GHz can be used either directly as dust tem- 
plates, or as extra information in a modified extrapolation 
of the dust to low frequencies, which will better trace the 
thermal dust seen at frequencies ~ 100 GHz. 



2.4 WMAP 7-year data 



We u se the l°-smoothed WMAP 7-year maps ( Jarosik et al 
201l[ ) provided by LAMBD^j^ which we then smooth fur- 
ther to the required resolution. This approach avoids prob- 
lems with e.g. non-circular beams, which will have been 
taken into account by WMAP. The maps are at frequen- 
cies of 22.8, 33.0, 40.7, 60.8 and 93.5 GHz. 

We also require the number of observations per pixel, 
Nohs, for the WMAP maps after they have been smoothed 
and their pixel resolution has been degraded. The effective 
A^obs maps were calculated using 1000 Monte-Carlo simu- 
lations. Each noise reaUzation, based on a Gaussian white 
noise with cr = ctq / \/Nohs where ao is the nominal thermal 
noise per observation (see Table 1), was smoothed to the 
appropriate resolution and degraded to the final A^side- The 
effective A'oba was then calculated using the r.m.s. scatter 
of the realizations for each pixel to get the final noise level. 



http : //lambda . gsf c. nasa . go v/product / map/dr4/ skymap_ 



and renormalizing to obtain the effective A^obs for each pixel 
for the maps as used in the analysis. 

We mask the Galactic plane using three different masks: 
the WMAP 7-year extended temperature mask, "KQ75", 
which masks 29.4 per cent of the sky; the WMAP 7-year 
temperature mask, "KQ85", which masks out 21.7 per cent 
of the sky ( Gold et al.|201l" I; and the "kp2" mask from the 
WMAP 1-year release, which masks only 15 per cent of the 
sky. Using these masks excludes areas of the sky that are 
very bright in emission that could bias the results, partic- 
ularly for fitting the free-free emission since Ha along the 
Galactic plane is absorbed by dust. In addition to mask- 
ing the galactic plane, all of these masks also exclude the 
pixels that contain bright extragalactic and other compact 
sources. Although there are some differences between the 1 
and 7 year point source masks, this will be small since these 
sources will be faint, and are further diluted by the map 
smoothing applied here. 



3 CROSS-CORRELATION METHOD 

We cross-correlate each WMAP map with the emission tem- 
plates by carrying out a least-squares fit to minimise the 
via a direct matrix inversion. 



as per 



Davies et al. (20061 



and Ghosh et al. (20111. The cross-correlation measure, a. 



linfo.cfml 



between a data vector, d and a template vector t can be 
measured by minimising: 
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X 



(d - at) 



(d - at) = d^ 



Ms^-d 



(1) 



where Msjv is the covariance matrix including both signal 
and noise for the template-corrected data vector d = d — at. 
Solving for a then becomes: 



■ Mg^ ■ d 

' ■ M.-ZI ■ t 



(2) 



To compare multiple template components tj, e.g. differ- 
ent foregrounds, to a given dataset, the problem becomes 
a matrix equation. In the case where we have different 
foreground components, we end up with the simple system 
of linear equations Ax = b, where 



Aki — tJ. ■ Mg-Jr ■ t, 



^SN ■ 

bk = t^ ■ Mc^T • d 



(3) 

(4) 
(5) 



When only one template is present, this reduces to equa- 
tion [2] above. 

We assess the uncertainty on our results using the full 
A^'pix X A^pix covariance matrijj^ The CMB anisotropies are 
taken into account statistically via the covariance matrix, 
calculated using the theoretical CMB anisotropy power spec- 
trumviaM,^ = ^YZoi'2^+'^)CtBjBl^,,Pi{n,-h,), where 
Be is the beam window function for a Gaussian beam of the 
appropriate resolution, and Be^pi^ is the HEALPix window 
function. The power spectrum d is the WMAP 7-year best- 
fit theoretical power spectrum ( Larson et al.|20TT l. 

We have verified that the uncertainty estimates from the 
cross-correlation method are robust, both when using the 
sync[H0.4] and sync[J2.3] templates, by using monte carlo 
simulations consisting of the foreground maps with realis- 
tic template coefficients (i.e. those given in Table [2]|; Gaus- 
sian noise based on WMAP noise levels and A^obs maps; and 
CMB realisations generated from the WMAP 7-year CMB 
power spectrum. Note that the dominant contribution to the 
covariance matrix is the CMB|3 



4 RESULTS 

Throughout this work, we only consider the area of the sky 
that has been surveyed at 2.3 GHz and lies outside of the 
Galactic plane mask; this is around half of the sky. We dis- 
cuss the basic results in Section |4.1[ followed by an inves- 
tigation of the robustness of these results against various 
parameters and options in Section [4. 2[ 

Throughout this paper, the units of the template co- 
efficients are \iK/K for sync[H0.4], [iK/mK for sync[J2.3], 
HK/R for Ha and [iK/[iK for dust[FDS94]. AU plots use 
a consistent colour scheme of black lines for sync, red 
lines for Ha and blue lines for dust, with solid lines for 



^ As A^pix can be up to ~30,000 for a masked map with A^sidc = 
64, this can require up to ~ 14 GB of RAM. 

^ Although these uncertainties could be significantly decreased 
by subtracting a map of the CMB from the data, this would 
introduce unquantified uncertainties and aliasing between com- 
ponents depending on the foreground residuals and noise in the 
CMB map. See Appendix A of |bavies et al., ( |2006| > for additional 
discussion of this issue. 
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Figure 3. Template coefficients for fits to all of the sky observed 
at 2.3 GHz outside of the KQ75 mask, using sync[H0.4] (solid 
lines) and sync[J2.3] (dashed lines). The red data points and 
lines are for Ha; blue arc for dust and black are for synchrotron. 
There is little difference between the solid and dashed lines 
for Ha[DDD] and dust[FDS94]. Note that only the top of the 
93.5 GHz error bar for the sync[H0.4] analysis is shown here as 
the central value is negative. 



sync[H0.4] and dashed lines for sync[J2.3]. For display pur- 
poses, we have renormalised the sync[J2.3] template coeffi- 
cients by (2.326/0. 408)~^ °/10"^ = 5.4, i.e. by assuming a 
spectral index of -3.0 between 408 MHz and 2.3 GHz. We 
have also offset the sync[J2.3] data points in frequency by 1 
per cent. All error bars are la. 



4.1 Basic results 

We apply the cross-correlation method to all 5 in- 
dividual WMAP frequency maps using the Ha[DDD], 
dust[FDS94] and sync[H0.4] templates, and also using the 
sync[J2.3] template in place of sync[H0.4]. The basic results 
for all of the sky observed at 2.3 GHz using maps smoothed 
to 3° and using A'side = 64, with the KQ75 mask applied, 
are shown in Figure [3j and are given numerically in Table [2] 
As expected, the synchrotron and free-free template co- 
efficients fall off as a function of frequency (i.e. they have 
negative spectral indices). The dust-correlated component 
decreases as a function of frequency where it is tracing the 



higher frequency side of the AME (see e.g. Planck Collab- 



oration et al. 20111; it then starts to increase at W-band 



where the thermal dust component begins to dominate. 

There are only a small changes in the Ha and dust 
correlation coefficients when the two different synchrotron 
templates are used, with the Ha coefficient systematically 
increasing by around 3 per cent between sync[H0.4] and 
sync[J2.3], and the dust coefficient systematically decreasing 
by around 7 per cent. This ~ 3cr decrease indicates that the 
synchrotron spectrum is flattening slightly at 2.3 GHz, how- 
ever it is not flattening sufficiently to explain all of the dust- 
correlated emission. For example, based on the rudimentary 
model described in the introduction, one would expect a ~50 
per cent decrease in order for the dust-correlated emission 
to be due to a flat spectrum synchrotron component with 
/3 = -2.5. 
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Table 2. Template coefficients using the sync[H0.4] and sync[J2.3] templates and the KQ75 mask (combined with the 2.3 GHz mask) 
for different regions: all-sky (top), North Galactic hemisphere (middle) and South Galactic hemisphere (bottom). The spectral indices 
are between the frequency in the row and the one before, except for the synchrotron spectral indices which are between frequency of the 
dataset and that of the synchrotron template. We omit the two highest frequency WMAP frequency bands as the synchrotron emission 
at those frequencies is low, and hence noisy. Note that the fractional uncertainty in the synchrotron template coefficients are consistent 
between the two synchrotron templates, even though the absolute values change. 



V 

(GHz) 


Ha 


sync[H0.4] 

Dust Sync /3sync 




/3dust 


Hq 


sync[J2.3] 

Dust Sync /3sync 


/3fi 


/^dust 


22.8 
33.0 
40.7 


10.0 ±0.3 7.9 ±0.2 6.7 ±0.2 
5.0 ±0.3 2.9 ±0.2 2.0 ±0.2 
3.2 ±0.3 1.7 ±0.2 0.9 ±0.2 


-2.96 ±0.01 
-2.99 ±0.03 
-3.02 ±0.05 


-1.9 ±0.2 
-2.0 ±0.5 


-2.7 ± 0.2 
-2.6 ±0.6 


10.3 ± 0.3 7.4 ± 0.2 0.83 ± 0.02 
5.1 ±0.3 2.7 ±0.2 0.29 ±0.02 
3.3 ±0.3 1.5 ±0.2 0.16 ±0.02 


-3.11 ± 0.01 
-3.07 ± 0.02 
-3.07 ± 0.04 


-1.9 ±0.2 
-2.0 ±0.5 


-2.7 ± 0.2 
-2.7 ± 0.6 


22.8 
33.0 
40.7 


9.4 ±0.6 

4.5 ±0.6 
2.8 ±0.6 


7.3 ±0.3 7.2 ±0.4 
2.7 ±0.3 2.3 ±0.4 
1.6 ±0.3 1.2 ±0.4 


-2.94 ± 0.01 
-2.95 ±0.04 
-2.97 ± 0.07 


-2.0 ±0.4 
-2.2 ± 1.1 


-2.7 ± 0.3 
-2.4 ± 1.1 


9.6 ±0.6 
4.6 ±0.6 
2.9 ±0.6 


6.6 ±0.3 0.93 ±0.03 
2.4 ± 0.3 0.35 ± 0.03 
1.4 ±0.3 0.20 ±0.03 


-3.06 ±0.01 
-3.00 ±0.03 
-2.98 ±0.04 


-2.0 ±0.4 
-2.2 ± 1.1 


-2.7 ± 0.4 
-2.5 ± 1.2 


22.8 
33.0 
40.7 


10.2 ±0.4 
5.1 ±0.4 
3.4 ±0.4 


8.0 ±0.2 6.1 ±0.4 
3.0 ±0.2 1.6 ±0.4 
1.7 ±0.2 0.5 ±0.4 


-2.98 ±0.02 
-3.04 ±0.06 
-3.15 ±0.16 


-1.8 ±0.2 
-2.0 ±0.6 


-2.7 ± 0.2 
-2.6 ±0.7 


10.1 ± 0.4 7.6 ± 0.2 0.76 ± 0.03 
5.1 ±0.4 2.8 ±0.2 0.26 ±0.03 
3.4 ±0.4 1.6 ±0.2 0.13 ±0.03 


-3.15 ±0.02 
-3.12 ±0.04 
-3.12 ±0.07 


-1.8 ±0.2 
-2.0 ±0.6 


-2.7 ±0.2 
-2.7 ±0.7 



The effect of the flat spectrum component at GHz fre- 
quencies can also be seen in the synchrotron spectral in- 
dices, where they are steeper between sync[J2.3] and the 
WMAP maps than they are with sync[H0.4]; this behaviour 
requires some fiat emission between 408 MHz and 2.3 GHz 
that falls off by WMAP frequencies. This is not too surpris- 
ing, since the average spectral index between sync[H0.4] and 
sync[J2.3] is /3 ~ -2.7 ( [Platania et al.||2003[ ). Given a 
spectral index between 408 MHz and 22.8 GHz of -2.96, 
a steeper index of ~ —3.16 would be expected between 
2.3 GHz and 22.8 GHz, i.e. similar to what we find here. Note 
that although calibration errors in the synchrotron maps will 
affect the spectral indices, any such errors would be expected 
to be around 10 per cent, and hence would not be an issue 
here. 

That the differences in the template coefficients when 
using the sync[J2.3] or sync[H0.4] templates are small con- 
firms that the bulk of the diffuse anomalous emission seen by 
WMAP is indeed dust-correlated, and hence is more likely 
to be spinning dust rather than synchrotron emission with 
a flattening spectral index. The dust-correlated emission at 
22.8 GHz is detected at a significance level of > 30a. 



4.2 Robustness 

To determine the robustness of our results, we look at the 
effects of different resolutions, foreground masks (including 
looking at the North and South Galactic hemispheres sepa- 
rately). Ha templates and using both synchrotron templates 
together. The results from these tests are described below, 
and are plotted in Figure |4] with the template coefficients 
presented in Table [XT] We find that our results as described 
above are robust. 

• Resolution: We have repeated the analysis above us- 
ing maps smoothed to 1°, 2° and 4° to look at the effects on 
the template coefficients due to different map resolutions; as 



shown by Ghosh et al. (2011 \ these can be considerable. The 
top-left panel of Figure |4| shows the template coefficients at 
22.8 GHz vs. the resolution of the map. Importantly, the dust 
coefficient remains relatively stable regardless of resolution. 



Table 3. Table of synchrotron spectral indices using the KQ75 
mask and resolutions of 1°, 2°, 3° and 4°. Top: using sync[H0.4]. 
Bottom: using sync[J2.3]. 



u (GHz) 


1° 


2° 


3° 


4° 


22.8 


-2.90 ±0.01 


-2.93 ±0.01 


-2.96 ±0.01 


-2.94 ± 0.01 


33.0 


-2.90 ±0.01 


-2.92 ±0.02 


-2.99 ±0.03 


-2.97 ±0.02 


40.7 


-2.91 ± 0.01 


-2.93 ±0.03 


-3.02 ±0.05 


-3.00 ±0.05 


22.8 


-2.85 ±0.01 


-2.87 ± 0.01 


-3.11 ± 0.01 


-3.13 ±0.01 


33.0 


-2.81 ± 0.01 


-2.84 ±0.01 


-3.07 ± 0.02 


-3.10 ±0.02 


40.7 


-2.79 ±0.01 


-2.82 ±0.01 


-3.07 ± 0.04 


-3.08 ±0.04 



The Ha template coefficient increases with larger beam size, 
with the Ha coefficients at 3° being consistent with a free- 
free electron temperature of 6000 — 7000 K, unlike previous 
studies such as'Davies et al. ( 2006 1. This is likely to be due to 



a number of effects including under-sampling reducing the 
smaller-scale power, stellar residuals, baseline offsets, etc.; 
see 



Ghosh et al. 



(2011 1 for discussion. 
The synchrotron coefficient decreases with larger beam 
size, potentially due to contaminating compact regions 
(which typically have flatter spectral indices) whose effect 
at larger resolutions has been reduced by the smoothing. 
The effect that this has on the synchrotron spectral index is 
shown in Table [3]- this can be quite large, particularly for 
2.3 GHz. There is a trend towards steeper spectral indices at 
coarser resolutions. These trends are present when using ei- 
ther sync[H0.4] or sync[J2.3]. The coefficients have however, 
converged by ~ 3° , verifying our choice of this resolution in 
the above analysis. 

• Mask: In order to check that our results are robust 
against the different masks that could be applied, and hence 
the different amount of strong Galactic plane emission that 
is masked out, we compare our results from the KQ75 mask 
at 22.8 GHz with those from the KQ85 and kp2 masks (as de- 
scribed in section [2!4| ). The results are shown in the top-right 
panel in Figure |4j The template coefficients increase as the 
size of the mask is decreased thereby including more of the 
Galactic plane. This is particularly the case for the Ha map. 




I 1 1 1 1 ' ' ' ' ' 

Ha[F03] Ha[DDD] Ha[DDD0.33] north both south 

Ha template Galactic hemisphere 

Figure 4. 22.8 GHz Template coefficients vs. resolution for all sky, showing convergence by 3° (top-left); vs. the WMAP galactic cut 
mask used (top-right); vs. the Ha template for all sky (bottom-left); and for all sky compared with North and South Galactic hemispheres 
(bottom-right). All panels use black lines for synchrotron, red lines for Ha and blue lines for dust, with solid lines indicating results using 
sync[H0.4], and dashed lines those using sync[J2.3]. The lines between data points are included to guide the eye as to the trends of the 
data rather than being physically meaningful. 



as would be expected if emission is being absorbed along the 
Galactic plane. There is a marked difference in the dust tem- 
plate for kp2 between sync[H0.4] and sync[J2.3], and a corre- 
sponding increase in the sync[J2.3] coefficient, which is likely 
due to the presence of flatter synchrotron spectra closer in 
to the Galactic plane. However, off the Galactic plane the 
results remain robust between sync[H0.4] and sync[J2.3]. 

We also consider the north- and south-Galactic hemi- 
spheres separately as well as combined; the results are shown 
in the bottom-right panel of Figure |4] and the numerical 
results are also given in Table |2] Although the results are 
largely consistent with each other, we find a slight asym- 
metry: the south Galactic hemispheres have a lower syn- 
chrotron correlation, and higher dust and Ha correlations, 
than the northern hemisphere. Although this is different 
from the results in [Ghosh et al.| ( |2011[ ), this is likely due 
to the exclusion of much of the Northern hemisphere from 
the analysis due to it not being covered by the 2.3 GHz map. 

• Ha template: We use the three different Ha templates 
to look for the dependency of our results on the accuracy 
of these maps. The results from analyses using both the 
HQ[F03]map and Ha[DDD] maps with and without dust 
corrections are shown in the bottom-left panel of Figure 



[4] We find that the dust correction significantly decreases 
the coefficient of free-free emission, and to a lesser extent 
the thermal dust emission, but that this equally affects the 
analyses using sync[H0.4] and sync[J2.3]. 



• Multiple synchrotron templates: We also use all 
four templates simultaneously to fit the WMAP data. Al- 
though the resulting synchrotron template coefficients are 
meaningless due to the large amount of degeneracy between 
the spatial structure in the two synchrotron templates, it 
is a useful cross-check for the Ha and dust coefficients. 
We find that this does not significantly affect our results. 
For example, at 22.8 GHz we find the template coefficient 
for Ha when using just sync[H0.4] to be 10.0 ± 0.3 ^.K/R 
and that for dust to be 7.9 ± 0.2 (i.K/[i.K. The coefficients 
change to 10.3 ± 0.3 ^.K/R and 7.3 ± 0.2 jiK/jiK respec- 
tively when including both sync[H0.4] and sync[J2.3], and 
to 10.3 ± 0.3M.K/R and 7.4 ± 0.2^K/^K when just using 
sync[J2.3]. In all cases, the dust coefficient varies by less 
than 10 per cent. 
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Figure 6. Template coefficients for fits using Haslam 408 MHz (solid lines) and Jonas 2.3 GHz (dashed lines). The red data points and 
lines are for Hq; blue are for dust and black are for synchrotron. Top-left is for the north polar spur; top-right for the Gum nebula; 
middle-left for Eridanus; middle-right for the Northern region; bottom for the South Galactic regions. 



4.3 Results in specific regions 

Although our results are robust on average across the en- 
tire sky, it is also important to look at their robustness in 
different Galactic environments. This also lets us test the 
effect of variations in the synchrotron spectral index across 
the sky. To do this, we focus on a set of five distinct regions, 
as outlined in Figure [S] in a similar way to |Davies et al 



(|2006[ ) . Since we are using a larger beam size than in 



Davies 



et al. (20061, we use significantly larger regions. This gives 



us more independent measures of the CMB fluctuations as 
well as allowing us to trace enough of the morphology to 
distinguish between the different components. [Ghosh et al.| 
(2011) look at a larger number of regions, based on the tem- 
plate morphology, which overlap with, but are distinct from, 
the regions we look at here. 

The results from applying template fitting to these five 
regions are shown in Figure [6j with the template coefficients 
presented in Table |A2| Although the uncertainties for each 
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Figure 5. Full-sky map showing the individual regions consid- 
ered at here, along with the 2.3 GHz survey area mask and the 
WMAP 7- year temperature mask (light dark). The regions are 
North Polar Spur (top; blue), Gum nebula (middle-right; red), 
Eridanus region (bottom-right and -left; yellow). South Galactic 
latitudes (bottom-middle; green) and the Northern region (top- 
right; orange) 

region are a lot higher, due to the smaUer area for each 
region, the dust and free-free coefHcients remain consistent 
within the uncertainties in all cases. The results from the 
individual regions are described below. 

• North Polar Spur (310° ^ I ^ 40°, 30° sC 6 sC 70°). 
This region is dominated by an arc of steep spectrum syn- 
chrotron radiation that is thought to be due to the remnant 
of a nearby supernova (see e.g. |Wolleben|2007| ) . The region 
also contains a spur of thermal dust emission (see Figure [2|. 

The results are shown in the top-left panel of Figure |6] 
Although there are large uncertainties on the template co- 
efficients, there is very little difference in the Ha and dust 
template coefficients between sync[EI0.4] and sync[J2.3]. 

• Gum nebula (240° ^ / 280°, -20° 6 ^ 20°). This 
region is dominated by free-free emission, with a reduced 
amount of synchrotron emission. It also contains a signifi- 
cant amount of thermal dust emission (see Figure |2| . We 
mask out the central part of the nebula, and look at the 
surrounding emission only, to avoid issues with the Ho? tem- 
plate in areas of high dust absorption. 

The results are shown in the top-right panel of Figure [5] 
The uncertainties are smaller since this is a bright region. 
We find a significantly increased synchrotron template co- 
efficient compared to the rest of the sky - approximately 4 
times higher - however this will be dominated by the un- 
certainty of the free-free correction to these templates. The 
Hq and dust coefficients between the two synchrotron tem- 
plates differ, with both being reduced when using sync[J2.3]. 
This implies that there is a reasonable amount of ffat spec- 
trum synchrotron emission in this region, but the result will 
depend on the free-free corrections to the synchrotron tem- 
plates. 

• Eridanus (170° < I 210°, -55° ^ 6 -25°). The 
main feature in this region is a line of thermal dust and free- 
free emission (as can be seen in Figure [2|, with a reduced 
amount of synchrotron emission compared to most of the 
sky 

The template coefficients are shown in the middle-left 
panel of Figure |6] The dust and Ha template coefficients 
from analyses using sync[H0.4] and sync[J2.3] are consistent 
with each other, however the uncertainties in this region are 
large as the foreground emission is small (in particular, the 



Table 4. Template cross-correlation matrices for the KQ75 mask 
(combined with the 2.3GHz mask), for all of the sky (top), the 
Northern region (middle) and the Eridanus region (bottom). 



Ha[DDD] dust[FDS94] sync[H0.4] sync[J2.3] const. 



i.±Lt LJ LJ 


1.00 


0.66 


64 


0.59 


67 


dust[FDS94] 


0.66 


1.00 


0.79 


0.75 


0.79 


sync[H0.4] 


0.64 


0.79 


1.00 


0.93 


0.96 


sync[J2.3] 


0.59 


0.75 


0.93 


1.00 


0.89 


const. 


0.67 


0.79 


0.96 


0.89 


1.00 


Ha[DDD] 


1.00 


0.40 


0.45 


0.39 


0.42 


dust[FDS94] 


0.40 


1.00 


0.61 


0.54 


0.63 


sync[H0.4] 


0.45 


0.61 


1.00 


0.86 


0.91 


sync [J 2. 3] 


0.39 


0.54 


0.86 


1.00 


0.80 


const. 


0.42 


0.63 


0.91 


0.80 


1.00 


Ha[DDD] 


1.00 


0.33 


0.26 


0.32 


0.28 


dust[FDS94] 


0.33 


1.00 


0.44 


0.47 


0.52 


sync[H0.4] 


0.26 


0.44 


1.00 


0.95 


0.88 


sync [J 2. 3] 


0.32 


0.47 


0.95 


1.00 


0.87 


const. 


0.28 


0.52 


0.88 


0.87 


1.00 



synchrotron emission is low) and hence the region is domi- 
nated more by CMB emission. 

• Northern region (220° ^ / 300°, 25° ^ fe 40°). 
This is a region of low foreground emission in the northern 
Galactic hemisphere, just above the Gum nebula. 

The results are shown in the middle-right panel of Figure 
[6] are consistent between sync[H0.4] and sync[J2.3](with a 
large uncertainty due to the low foregrounds), with a slightly 
higher dust coefficient when sync[J2.3] is used rather than 
sync[H0.4]. 

• South Galactic latitudes (320° ^ Z s: 70°, -60° s: 
b ^ —20°). This is a region to the south of the Galactic 
centre that contains moderate levels of foreground emission, 
primarily synchrotron and thermal dust emission, with a 
lower level of free-free emission. Figure [2] shows that this 
region has not been well-mapped by Ha surveys, with clear 
mottling present due to the field-of-view of the observations. 

The results are shown in the bottom panel of Figure [6] 
We find coefficients for the Ha template that are consistent 
with zero within large uncertainties, confirming the presence 
of issues with the Ha template in this region. There are high 
dust and synchrotron coefficients, which are consistent with 
each other within the uncertainties. 

4.4 Template cross-correlations 

Correlations between templates can result in degeneracies 
between the different template coefficients, which are im- 
portant to consider in this type of analysis. The correlations 
can be calculated from the off-diagonal terms of A from 
Equation [3] via 

Cij — Aij / y^AnAjj. (6) 

Table |4] gives the correlation matrices between the different 
types of templates for the whole of the area of sky used 
here, as well as in two illustrative regions (Northern and 
Eridanus) . 

As expected, in all cases the two synchrotron templates 
are highly correlated with each other, but not completely. 



10 M. W. Peel et al. 



Typical values for the correlation between the synchrotron 
templates are around 0.9. On the whole sky, there are sig- 
nificant correlations between all of the templates, due to 
the similar large-scale morphology of the Galactic plane and 
lower latitude regions. There is also a high degree of correla- 
tion with the constant template. The correlations decrease 
when looking at smaller regions of sky, where the large-scale 
structure of our Galaxy is less significant. That the cor- 
relation coefficients vary between different regions, whilst 
the regions have the same results, indicates that these cor- 
relations do not present a significant issue. This issue has 
been explored in more detail in Appendix B of |Davies et al.| 
(20061). 



An improved approach to take here would be to filter 
the maps and the templates to remove large-scale modes, 
and to concentrate on the smaller scale fluctuations, which 
will be investigated in a future work. 



5 CONCLUSIONS 

We have applied the template fitting method to the 
WMAP 7-year maps, using synchrotron templates at 
408 MHz and 2.3 GHz. We find that using the higher fre- 
quency synchrotron template slightly reduces the dust- 
correlated component across the whole sky by up to 10 
per cent, which indicates that the bulk of the emission is 
much more likely to be due to spinning dust emission rather 
than fiattening synchrotron radiation, in agreement with 



the results of e.g. de Oliveira-Costa et al. (20041; Planck 



Collaboration et al. (2011 1. We find that our results are ro- 



bust against different choices of the map resolution. Galactic 
masks and Ha template, and that our uncertainty estimates 
are robust. We find a slight asymmetry in the template coef- 
ficients between the north and south Galactic hemispheres, 
and also different template coefficients in different regions of 
the sky. However we find that the dust coefficient remains 
robust when using either synchrotron template in all cases 
off the Galactic plane. 

In the future, data from Planck LFI will provide addi- 
tional measurements that can be compared with the tem- 
plates to help pin down the spectra of the different com- 
ponents. A major step forward will be when the 5 GHz 
all-sky map currently being produced by C-BASS becomes 
available; as this is all-sky, and at twice the frequency of 
sync[J2.3], this will greatly improve the constraints on flat 
synchrotron emission at higher frequencies than those con- 
sidered here, and will further help pin down the cause of the 
anomalous microwave emission. 
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sync[H0.4] sync[J2.3] 



Identifier 


Ho 


Dust 


Sync 


Z^sync 


Ha 


Dust 


Sync 


/^sync 


1° 


6.2 ±0.1 


7.6 ±0.2 


8.5 ±0.1 


-2.90 ±0.01 


6.9 ±0.1 


7.3 ±0.2 


1.5 ± 0.1 


-2.85 ±0.01 


2° 


9.2 ± 0.2 


7.8 ± 0.3 


7.7 ±0.1 




9.7 ± 0.3 


7.0 ± 0.3 


1.4 ± 0.1 


ij.Ol _i_ W.W-L 


3° 


10.0 ±0.6 


7.9 ±0.3 


6.7 ± 0.2 


-2.96 ±0.01 


10.3 ±0.7 


7.4 ± 0.3 


0.8 ±0.2 


-3.11 ± 0.01 


4° 


10.6 ±0.8 


8.1 ± 0.3 


7.3 ±0.2 


-2.94 ± 0.01 


10.7 ± 0.8 


7.8 ±0.3 


0.8 ±0.2 


-3.13 ±0.01 


KQ75 


10.0 ±0.6 


7.9 ±0.3 


6.7 ± 0.2 


-2.96 ±0.01 


10.3 ±0.7 


7.4 ± 0.3 


0.8 ±0.2 


-3.11 ± 0.01 


KQ85 


10.4 ± 0.9 


8.3 ±0.2 


7.3 ±0.1 


-2.94 ± 0.01 


10.4 ±0.9 


8.2 ±0.2 


0.8 ±0.1 


-3.11 ± 0.01 


kp2 


12.0 ±0.8 


12.6 ±0.1 


6.5 ±0.1 


-2.97 ± 0.01 


12.1 ± 0.8 


10.1 ± 0.1 


1.8 ±0.1 


-2.78 ±0.01 


Ha[F03] 


10.1 ± 0.6 


7.8 ±0.3 


6.7 ± 0.2 


-2.96 ±0.01 


10.3 ±0.7 


7.3 ±0.3 


0.8 ±0.2 


-3.11 ± 0.01 


Ha[DDD] 


10.0 ±0.6 


7.9 ±0.3 


6.7 ± 0.2 


-2.96 ±0.01 


10.3 ±0.7 


7.4 ±0.3 


0.8 ±0.2 


-3.11 ± 0.01 


Ho[DDD0.33] 


8.4 ± 0.6 


7.1 ± 0.3 


6.7 ± 0.2 


-2.96 ±0.01 


8.7 ± 0.7 


6.7 ± 0.3 


0.8 ±0.2 


-3.11 ± 0.01 


north 


9.4 ± 0.4 


7.3 ±0.6 


7.2 ±0.3 


-2.94 ± 0.01 


9.6 ±0.5 


6.6 ±0.6 


0.9 ±0.3 


-3.06 ±0.01 


both 


10.0 ±0.6 


7.9 ±0.3 


6.7 ± 0.2 


-2.96 ±0.01 


10.3 ±0.7 


7.4 ±0.3 


0.8 ±0.2 


-3.11 ± 0.01 


south 


10.2 ±0.5 


8.0 ±0.4 


6.1 ± 0.2 


-2.98 ±0.02 


10.1 ± 0.6 


7.6 ±0.4 


0.8 ±0.2 


-3.14 ± 0.02 



Table A2. Template coefficients, as per Table [2] for the five regions shown in Figure |6] From top to bottom: North Polar Spur, Gum 
nebula, Eridanus, Northern and South Galactic. 



V 

(GHz) 


Ha 


Dust 


sync[H0.4] 
Sync ftync 




/^dust 


Ha 


Dust 


sync[J2.3] 
Sync Aync 




/^dust 


22.8 
33.0 
40.7 


13.6 ± 2.2 
8.1 ± 2.1 
6.3 ±2.1 


9.5 ± 0.8 
3.0 ±0.8 
1.4 ±0.8 


5.0 ±0.7 
1.2 ±0.7 
0.3 ±0.7 


-3.03 ±0.03 
-3.10 ±0.13 
-3.25 ±0.47 


-1.4 ±0.8 
-1.2 ± 2.0 


-3.1 ± 0.7 
-3.5 ± 2.8 


11.4 ±2.2 
7.2 ± 2.2 
5.6 ± 2.1 


9.2 ±0.8 
2.9 ±0.8 
1.4 ± 0.8 


0.81 ± 0.09 
0.33 ±0.09 
0.20 ± 0.09 


-3.12 ±0.05 
-3.02 ±0.10 
-2.97 ± 0.15 


-1.3 ± 1.0 
-1.1 ± 2.3 


-3.1 ± 0.8 
-3.5 ±2.9 


22.8 
33.0 
40.7 


14.6 ±0.1 
6.9 ±0.1 
4.3 ±0.1 


12.8 ±0.1 
5.7 ±0.1 
3.5 ± 0.1 


52.9 ±0.3 
18.4 ± 0.3 
9.9 ±0.3 


-2.45 ± 0.00 
-2.48 ±0.00 
-2.50 ±0.01 


-2.0 ±0.0 
-2.2 ± 0.1 


-2.2 ± 0.0 
-2.3 ±0.1 


13.1 ± 0.1 
6.2 ±0.1 
4.0 ±0.1 


10.4 ± 0.1 
4.5 ±0.1 
2.7 ± 0.1 


4.17 ± 0.01 
1.72 ±0.01 
1.03 ±0.01 


-2.40 ±0.00 
-2.40 ±0.00 
-2.40 ±0.00 


-2.0 ±0.1 
-2.1 ± 0.2 


-2.3 ±0.0 
-2.4 ±0.1 


22.8 
33.0 
40.7 


8.2 ± 0.5 
4.4 ±0.5 

3.3 ±0.5 


9.8 ±0.4 
3.7 ±0.4 

1.9 ±0.4 


9.2 ±2.4 
4.8 ± 2.4 
4.0 ± 2.4 


-2.88 ±0.07 
-2.79 ±0.11 
-2.70 ±0.13 


-1.7±0.4 
-1.4± 1.0 


-2.6 ±0.3 
-3.0 ± 1.2 


8.0 ±0.5 
4.3 ±0.5 

3.1 ± 0.5 


10.1 ± 0.4 
3.8 ±0.4 
2.1 ± 0.4 


0.58 ±0.12 
0.35 ±0.11 
0.25 ±0.11 


-3.27 ± 0.09 
-3.00 ±0.12 
-2.90 ±0.16 


-1.7 ±0.4 
-1.4± 1.0 


-2.6 ±0.3 
-2.9 ±1.0 


22.8 
33.0 
40.7 


12.0± 3.1 
7.9 ±3.0 
6.2 ± 3.0 


5.8 ± 1.9 
0.4 ± 1.9 
-0.7± 1.8 


10.2 ±2.0 
6.2 ± 1.9 
5.2 ± 1.9 


-2.86 ±0.05 
-2.73 ±0.07 
-2.64 ±0.08 


-1.1 ± 1.3 
-1.1 ±2.9 


-7.1 ± 12.0 


13.0 ±3.1 
8.9 ±3.0 
7.3 ±3.0 


6.1 ± 1.9 
0.7 ± 1.9 
-0.3 ± 1.8 


0.94 ± 0.20 
0.32 ±0.19 
0.15 ±0.19 


-3.05 ±0.09 
-3.03 ±0.23 
-3.09 ± 0.45 


-1.0± 1.1 
-1.0 ±2.5 


-5.7 ± 6.9 


22.8 
33.0 
40.7 


1.5 ± 2.9 
1.0 ±2.9 
-0.1 ± 2.9 


7.7 ±0.7 
2.0 ±0.7 
0.7 ±0.6 


6.8 ±0.9 
1.8 ±0.8 
0.6 ±0.8 


-2.96 ±0.03 
-3.02 ± 0.11 
-3.13 ±0.32 


-1.3±9.7 


-3.6 ±0.9 
-5.1 ±4.7 


2.2 ± 2.9 
1.0 ± 2.9 
-0.2 ±2.8 


7.7 ± 0.7 
2.0 ± 0.7 
0.6 ±0.6 


0.97 ± 0.08 
0.32 ±0.08 
0.16 ±0.08 


-3.04 ±0.04 
-3.03 ±0.10 
-3.06 ±0.18 


-2.3 ±8.9 


-3.7 ±0.9 
-5.4 ±5.1 
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